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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase, a superconducting state with non-zero total 
momentum Cooper pairs in a large magnetic field, was first predicted about 50 years ago, and since 
then became an important concept in many branches of physics. In recent years, the possibility 
of observing FFLO states using ultracold degenerate Fermi gases has sparked tremendous interest. 
However, unambiguous experimental evidence for FFLO states is still elusive because of the stringent 
parameter requirement in experiments. In this Letter, we show that a giant parameter region for 
FFLO states can be obtained in 3D degenerate Fermi gases in the presence of spin-orbit coupling 
and an in-plane Zeeman field, two ingredients that were already developed for cold atoms in recent 
experiments. The predicted FFLO state is stable against quantum fluctuations due to the 3D 
geometry, and can be observed with experimentally already achieved temperature (T ~ 0.05-Ef), 
thus opens a new fascinating avenue for exploring FFLO physics in degenerate Fermi gases. 
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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase, 
characterized by Cooper pairs with nonzero total momen- 
tum and spatially non- uniform order parameters, was 
predicted to exist in certain region of superconductors 
in high Zeeman fields [2-01 . This fascinating state arises 
from the interplay between magnetic and superconduct- 
ing order, and now is a central concept for understand- 
ing many exotic phenomena in different physics branches, 
ranging from unconventional solid state superconductors 
(e.g., layered 0, 0], heavy-fermion 0-0] ; organic (il Hoj 
superconductors, etc.), to chiral quark matter in quan- 
tum chromodynamics , and to neutron star glitches in 
astrophysics [ill, E3] ■ Despite tremendous experimental 
and theoretical efforts in the past five decades, there is 
still no unambiguous experimental evidence for FFLO 
states [12] • The experimental difficulty may arise from 
several different aspects, such as the depairing of Cooper 
pairs due to orbital or Pauli paramagnetic effects in 
strong magnetic fields and unavoidable disorder effects 
in solid state materials. 

The recent experimental realization of spin-imbalanced 
Fermi gases [1314171 ] provides a new excellent platform for 
exploring FFLO physics. In Fermi gases, the effective 
Zeeman field is generated through the population imbal- 
ance between two spins, therefore the orbital effects (e.g., 
vortices induced by the magnetic field) are absent even 
in 3D. The Fermi gases are also free of disorder and all 
experimental parameters are highly controllable. These 
advantages have sparked tremendous recent interest in 
exploring FFLO physics in spin-imbalanced Fermi gases 
[18l426| |. However, the FFLO phase only exists in a nar- 
row parameter regime in 3D due to the Pauli param- 
agnetic depairing effect 18, 22, |2j|. Furthermore, the 



transition from the BCS superfluid to the normal gas 



energy difference between the FFLO states and the BCS 
superfluid is extremely small [HI- As a result, only the 



nun I has been observed in experiments in 3D spin- 
imbalanced Fermi gases. Current experimental and the- 
oretical efforts on the FFLO state have focused on low 
dimensions (ID or 2D) [27T - [3lj | . where quantum and ther- 
mal (at finite temperature) fluctuations may become cru- 
cial and the physics is much more complicated. 

In this Letter we show that a large and stable param- 
eter region for FFLO states can be realized even in 3D 
degenerate Fermi gase s by including two experimentally 
already developed [32j-|34| elements: spin-orbit (SO) cou- 
pling and an in-plane Zeeman field. Recently the prop- 
erties of SO coupled Fermi gases with perpendicular Zee- 
man fields have been intensively investigated with the 
goal of realizing topological superfluids [35l - l38j | and the 
associated Majorana fermions 3^-4l|. However, regular 
BCS superfluids, instead of FFLO states, are energeti- 
cally preferred in the presence of perpendicular Zeeman 
fields because of the centrally symmetric Fermi surface. 
We show that this issue can be resolved by using an in- 
plane Zeeman field, which, together with the SO cou- 
pling, yields an asymmetric Fermi surface so that the 
FFLO state can emerge naturally. At the same time, the 
Pauli paramagnetic depairing effect can be greatly sup- 
pressed by the SO coupling, leading to a large parameter 
region for the FFLO state. More importantly, we find 
that the energy difference between the FFLO ground 
state and the possible BCS superfluid excited state is 
dramatically increased (to ~ O.OAEp per particle), there- 
fore the FFLO state is experimentally more accessible 
with the realistic temperature in 3D (T ~ Q.05Ef). Fi- 
nally, because of the 3D geometry, the quantum and ther- 
mal fluctuations that play major roles in ID and 2D arc 
strongly suppressed, which greatly simplifies the FFLO 
physics. 



2 




-1.0 -0.5 0.0 0.5 1.0 



1/K F a s 



-1.0 -0.5 0.0 0.5 1.0 



1/K F a s 




-1.0 -0.5 0.0 0.5 1.0 



1/K F a s 



ocKp/Ep 



FIG. 1: Phase diagrams in the presence of SO coupling and 
in-plane Zeeman field, (a) Without SO coupling. The circle 
symbol represents the data from the quantum Monte Carlo 
calculation [4Q]. (b) and (c) With SO coupling olKf — O.HEf 
and olKf = I.OEf- (d) In the unitary regime (1/Kfcls = 
0.0). 



Consider a 3D degenerate Fermi gas in the presence of 
a Rashba type of SO coupling and an in-plane Zeeman 
field. The corresponding partition function of the sys- 
tem can be expressed as Z = Tre~' 3 (^ f ~' lAr ) = J Vipe~ s , 
with the action S = J ip^(d T + H )i/j + gi/iji/'jV'iV't- Here 
/ = r,f drd 3 r, = H = £L-f,-ha x + 

X Gy P y <T X ) , m is the mass of the atom, [i is the chem- 
ical potential, g is the s-wave interaction strength, a is 
the Rashba SO coupling strength, and h is the in-plane 
(same as the SO coupling) Zeeman field. In experiments, 
the SO coupling and the in-plane Zeeman field can be re- 
alized using the tripod scheme where three Raman lasers 
couple three hyperfine ground states with a common ex- 



cited state |42h44| . Note that without SO coupling, the 
Zeeman field can be applied through adjusting the pop- 
ulation imbalance between two spin states. While in the 
presence of SO coupling, the number of atoms at each 
spin channel is not conserved anymore, and the popula- 
tion imbalance can be tuned only through the external 
Raman lasers (Zeeman field). In addition, an in-plane 
Zeeman field is generated naturally using three Raman 
lasers in the tripod scheme (42l - |44| . while a perpendicu- 
lar Zeeman field requires additional lasers [i^ (thus more 
difficult in experiments). 

In the FFLO state the Cooper pairs have finite total 
momentum, i.e., A(r) = (V-'iV't) = Ae iCJr , where Q is 
the FFLO vector. We adopt a spatial uniform order pa- 
rameter A in our calculation through a transformation 
of the field ip ~~ > ipe tC *' r / 2 , yielding a new Hamiltonian 



e iQ-r/2# ( p ) e iQ-r/2 = H ^ p + q/ 2 ) = R ^ and the &Q _ 

tion 

3= I V t (3 T + 5o)^-|A| 2 / 3 + A^V| + A t ^ t (!) 



in the mean field approximation. Integrating out the 
Fermi field, we obtain Z = J T>A exp(— S c s), with the 
effective action 
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where^q_ kff = (§-k) 2 /2m-^, E\ (A = I, 2, 3, 4) are 
the eigenstates of the matrix (under the basis (^Q/2+p.t! 
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In Eq. © the bare interaction strength g should be 
regularized in terms of the s-wave scattering length a s 
51 EH, i — \ — = 1 + y\ T~, where e k = 

The ground state phase diagram of the system (i.e., A, 
/j, Q) is determined by the saddle point of the thermo- 
dynamical potetntial ^ = and ^ = 0, as well as the 

atom number conservation n = X) CT =t I na = ~^7' wnere 
il = S c h//3. The energy unit is chosen as the Fermi en- 
ergy Ep for an non-interacting gas without SO coupling 
and Zeeman field. The length unit is Kp 1 . We restrict to 
T = throughout this work. Generally the vector Q has 
three different components, and the total five unknown 
parameters put a great burden for numerically solving the 
above equations self-consistently because the landscape 
of f2 is an extremely complex function of these parame- 
ters whose global minimum (instead of a local minimum) 
is hard to find. For the x-axis Zeeman field and the 
Rashba-type SO coupling the deformation of the Fermi 
surface is along the y-axis, therefore the FFLO vector is 
expected to be along the y axis, i.e., Q = (0, Q,0). We 
have numerically confirmed that there is no large FFLO 
region when Q is along the x and z directions. There 
are three possible phases in this system: BCS supcrfluid 
(A ^ 0, Q = 0), FFLO (A ^ 0, Q ^ 0), and normal gas 
(A = and Q = 0). In the FFLO phase, we also cal- 
culate the energy difference between the FFLO ground 
state and the possible BCS superfluid excited state (by 
forcing Q = 0) to check the stability of the FFLO state 
against the finite temperature effect. 

In Fig. [TJ we plot the phase diagrams of the Fermi gas 
with respect to the Zeeman field h, the s-wave interaction 
1/Kpds, and the SO coupling strength aKp. Without 
SO coupling (Fig. QJi), our numerical result agrees well 
with that in pre vious literature using the mean field ap- 
proximation [lc| or quantum Monte Carlo (46[. We see 
the FFLO phase exists within an extremely small regime 
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FIG. 2: Illustration of the physical mechanism of the FFLO 
state in the presence of an in-plane Zeeman field and SO cou- 
pling. Solid and dashed contours are two Fermi surfaces. The 
solid arrows are the pseudospins. The solid line connecting 
two pseudospins represents the Cooper pair with total mo- 
mentum Q (direction is shown by the dashed arrow), (a) 
Without SO coupling, the Fermi surfaces are two concen- 
tric spheres, (b) With SO coupling, the Fermi surfaces are 
anisotropic along the k y axis due to the Rashba SO coupling 
and the a;-axis Zeeman field. 



in the phase diagram without the SO coupling. Fur- 
thermore, the energy difference per particle between the 
FFLO state and the possible BCS superfluid state (ob- 
tained by forcing Q = 0) is extremely small (see Fig. [3jl), 
therefore the Fermi gas may not relax to the FFLO state 
considering the realistic temperature in experiments [3] , 
even if the FFLO state is the true ground state. With in- 
creasing SO strength, the parameter region of the FFLO 
phase is greatly enlarged. In Fig. [TJl, we see the the main 
enlargement of the FFLO phase comes from the decrease 
of the critical Zeeman field for the transition from the 
BCS superfluid to the FFLO phase. 

The enlarged parameter region for the FFLO state in 
Fig. Q] can be understood from the change of the shape 
of the Fermi surface due to the SO coupling and the 
in-plane Zeeman field. Without the SO coupling, the 
Zeeman field (no matter which direction) yields two con- 
centric spheres (Fig. [5^) of the Fermi surface, and only 
singlet pairing between different pseudospins (i.e., two 
eigenstates of H ) is allowed due to the SU(2) symmetry 
of the Hamiltonian. With increasing Zeeman fields, the 
Fermi surface mismatch increases the energy cost of the 
BCS superfluid. In a strong Zeeman field the superfluid 
has to break the spatial symmetry to lower the accumu- 
lated energy, therefore the FFLO state emerges, but only 
in a small parameter region due to the Pauli paramag- 
netic depairing effect. Such depairing effect in strong Zee- 
man fields can be circumvented using the SO coupling, 
which allows both singlet and triplet pairings (361 1471.148 1 
(the later is insensitive to the depairing effect) because 
the pseudospin state is a sp in mixed state with strong 
momentum dependence [35|. However, if a perpendicu- 
lar Zeeman field is applied, the regular BCS superfluid in 
the same SO band (i.e., triplet pairings) dominates be- 
cause of the symmetric Fermi surface, and the parameter 
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FIG. 3: BEC-BCS crossover in the presence of SO coupling 
and an in-plane Zeeman field, h — 0.5-Ef and oKf = 0.0 
and 1.0-Ep. In (a) and (b), the solid lines are obtained by 
minimizing the total free energy with respect to A, ji and Q, 
while the dashed lines are obtained by forcing Q = (thus no 
FFLO states), (c) Plot of Q as a function of the scattering 
interaction, (d) The free energy difference between the FFLO 
state and the possible BCS superfluid. Here the free energy 
F = Q + fin. 



region for the FFLO state indeed shrinks comparing with 
that with only Zeeman fields, as found in our numerical 
simulation. In contrast, in the presence of SO coupling 
and an in-plane Zeeman field, the Fermi surfaces become 
anisotropic and the center of the Fermi surface is also 
shifted accordingly (Fig. [5}d). Therefore the regular BCS 
superfluid, which is preferred for a symmetric Fermi sur- 
face, is greatly suppressed, and the FFLO state becomes 
energetically favorable in a much wider parameter region, 
as observed in Fig. [T] 

To characterize the FFLO state, in Figs. [3^ and[3}}) we 
plot the chemical potential \x and the order parameter A 
in the BCS-BEC crossover. For comparison, we also plot 
/i and A for the possible BCS superfluid state (by forcing 
Q = 0). In the weak BCS limit A is exponentially small, 
therefore a small Zeeman field or population imbalance 
can destroy the superfluid 



14J. In the BEC side, the 



fcrmions form tightly bound molecules and the influence 
of Zeeman field and SO coupling is negligible. There- 
fore the only relevant parameter regime for the observa- 
tion of FFLO states should be near the unitary regime. 
In the FFLO regime, A for the FFLO state is smaller 
than that for the assumed BCS superfluid to reduce the 
FFLO energy. In Fig. [3b, we plot Q versus the scatter- 
ing interaction, which also confirms that the SO coupling 
can greatly increase the parameter region for the FFLO 
phase. 
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FIG. 4: Plot of the population imbalance P = Sn/n as a 
function of the scattering interaction (a) and the SO coupling 
strength (b). 



An experimentally observable FFLO state requires a 
large energy difference between the FFLO ground state 
and the possible BCS superfluid excited state so that the 
FFLO state can survive at finite temperature. In Fig. 
[3JJ, we plot the difference between the free energies of the 
FFLO state and the BCS superfluid per particle, 8F — 
(-Ffflo — -Fbcs) /nE F , with F = Q + [in. The stability 
of the FFLO state has not been emphasized in previous 
literatures [H [H-El Hf| . For FFLO states without SO 
coupling we find SF ~ 10~ 4 Ep, which is much smaller 
than the experimental temperature (T ~ 0.05Ep) [13ll49| | 
by more than two orders of magnitude. Therefore the 
FFLO state cannot be observed even the exact parameter 
region has been reached. While with the SO coupling and 
in-plane Zeeman field, the energy difference is greatly 
enhanced to ~ 0.04Ep, which makes the FFLO state 
accessible with realistic experimental temperature. Such 
a large energy difference is another major advantage of 
our scheme over previous Zeeman field 0S] or optical 
lattice [2J, [26J schemes. 

The FFLO state has a strong relation with the pop- 
ulation imbalance of the Fermi gas. Since the Zeeman 
field is applied along the x axis, we define the population 
imbalance as P — 5n/n with Sn — (a x ). In Fig. Hk, 
we plot P with respect to 1/Kpa s . Without SO cou- 
pling, the BCS superfluid breaks down at P ~ 0.669, in 
consistent with previous results [l3|, [l8| . When the SO 
coupling is applied, the transition from the FFLO state 
to the normal gas is expected to occur at a larger pop- 
ulation imbalance due to the induced triplet pairings in 
the FFLO states which are insensitive to the depairing 
effects. In Fig. |4Jd, we plot P with respect to aKp in 
the unitary region. We see the population imbalance is 
greatly enhanced by the SO coupling because the Zee- 
man field and the SO coupling lie in the same plane. 
The same reason also leads to the reduced critical Zee- 
man field for the transition between the FFLO state and 
the BCS superfluid in Fig. [TJi. 

So far we only consider the existence of the FFLO state 
using a simplified pairing model A(x) — Ae iC ^' r , while 
the true pairing of the FFLO state may be different. Be- 



cause the FFLO state depends strongly on the nesting of 
the Fermi surface, the order parameter may be composed 
of multiple vectors @, Hoj], i.e., A(r) = A(Qi, Q 2 , ), 
whose stability depends strongly on the detailed struc- 
ture of the Fermi surface and thus cannot been ruled out 
[50| . However, the main conclusion of our work, the large 
parameter region and the stable FFLO state induced by 
SO coupling and in-plane Zeeman field, is in intact even 
for very complex pairings because different choices of the 
order parameter are mainly used to further reduce the to- 
tal energy of the FFLO state (thus further enhance our 
results). 

The FFLO states may be measured using the time- 
of- flight imaging [24l|. where the momentum distribution 
is symmetric for the BCS superfluid, but asymmetric for 
the FFLO state due to the finite Q. The spatial variation 
of the order parameter (thus the atom density) may also 
be observed using the phase-contrast image. The super- 
fluidity of the FFLO states can be demonstrated through 
the rotation of the system [5l| , where the generated vor- 
tices provides the signatures of superfluidity. Note that 
near the boundary of different phases the vortices may 
be unstable due to strong damping effects however 
in the middle of the FFLO phase (only possible with a 
large parameter region for the FFLO state), we expect 
the damping effect to be small, similar as that in the BCS 
superfluid state. 

Finally, we briefly comment on the FFLO state in 2D 
Fermi gases. In solid state systems, 2D superconducting 
materials are generally used for the observation of the 
FFLO state due to the absence of orbital effects when an 
in-plane magnetic field is applied. However, in most su- 
perconducting materials, the SO coupling is absent, and 
the observation of the FFLO state is thus difficult. We 
study the FFLO phase in 2D Fermi gases with the same 
model Hamiltonian at zero temperature, and find that 
the FFLO regime is also greatly enhanced. However, at 
finite temperature, the phase fluctuation effect is signif- 
icant, therefore the long-range order parameters do not 
exist and the relevant physics is the Kosterlitz-Thouless 
transition [39]. 

In summary, we show that the combination of SO 
coupling and in-plane Zeeman field can lead to a large 
and stable parameter region for the experimentally long- 
sought FFLO state even for 3D degenerate Fermi gases. 
Considering the recent experimental progress on the gen- 
eration of the SO coupling in Bose and Fermi gases, our 
work provides a new exciting research direction for the 
study of SO coupled Fermi gases as well as the FFLO 
physics, which is essential for the understanding of im- 
portant phenomena in many branches of physics, ranging 
from solid state superconductors to astrophysics. 
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